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Abstract 

A new method for generating exactly solvable Schròdinger equations with a 
position-dependent mass is proposed. It is based on a relation with some deformed 
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Schròdinger equations, which can be dealt with by using a supersymmetric quantum 
mechanical approach combined with a deformed shape-invariance condition. The solv- 
ability of the latter is shown to impose the form of both the deformed superpotential 
and the position-dependent mass. The conditions for the existence of bound states 
are determined. A lot of examples are provided and the corresponding bound-state 
spectrum and wavefunctions are reviewed 

1 Introduction 

Schròdinger equations with a position-dependent mass (PDM) play an important role in 
many physical problems. They appear in the energy-dependent functional approach to 
quantum many-body systems [1] (e.g., nuclei, quantum liquids, 3 He clusters, metal clus- 
ters) and are very useful in the description of electrònic properties of condensed-matter sys- 
tems [2] (e.g., compositionally-graded crystals, quantum dots, liquid crystals). The PDM 
presence may also reflect other unconventional effects, such as deformation of the canonical 
commutation relations or curvature of the underlying space [3] or else pseudo-Hermiticity 
of the Hamiltonian [4] . 

Several exactly solvable, quasi-exactly solvable or conditionally exactly solvable PDM 
Schròdinger equations have been constructed using point canonical transformations, Lie al- 
gebraic methods or supersymmetric quantum mechanical (SUSYQM) and shape-invariance 
(SI) techniques (see [5, 6] and references quoted therein). Most of them can be obtained 



1 



from known constant-mass models by changes of variable and of function. As a conse- 
quence the spectrum is left unchanged although the potential is given by a complicated 
mass-deformed expression. 

Here we present a new method [7], generalizing the work in [3]. It is based on a re- 
lation with deformed Schròdinger equations and uses a SUSYQM method combined with 
a deformed SI condition. The resulting spectra will contrast with those obtained in the 
constant-mass case. 



2 PDM and deformed Schròdinger equations 



As is well known, the noncommutativity of a PDM with the momentum operator results in 
an ambiguity in the kinetic energy operator definition coming from the selected ordering. On 
choosing the von Roos general two-parameter form of T [8] , which has an inbuilt Hermiticity 
and contains all the plausible forms as special cases, and setting h = 2m = 1, the PDM 
Schròdinger equation can be written as 



1 ( m ?±m*±a{C + Af? ±m* ± M e) + V 

2 \ dx dx dx dx ) 



Í> = Eip. 



Here V(o; x) is the potential, M(a; x) is the dimensionless form of the mass function 
m(ct; x) = m M(o:; x), a and ex denote two sets of parameters, and the von Roos ambiguity 
parameters ?/, Ç are constrained by the condition £' + rf + Ç — — 1. 

On setting M(ct;x) = [f (ct; x)]~ 2 and f(ct;x) = 1 + g(a;x), where f(ct;x) is some 
positive-definite function and g(oc;x) = corresponds to the constant-mass case, Eq. (1) 
becomes 

4 (f^r^f+f^r^A+v 



dx dx' 



dx dx" 



(2) 



with £ + r] + ( = 2. The ambiguity parameters £, r], ( (denoted collectively by £) can be 
transferred from the kinetic energy term to an effective potential 



V cS (b;x) = V + V, 
so that Eq. (2) aequires the form 



V = pff" + of\ 



^4 



+ v c 



cS 



(3) 



(4) 



In (3), a prime denotes derivative with respect to x, p — |(1 — £ — (), a = (| — £) (| — C) 
and b = (a, oc, £). 

Equation (4) can be reinterpreted as a deformed Schròdinger equation with H = 7T 2 +V e e, 
where ir = \ff V^fl satisfies the deformed commutation relation [x, n] = if and /(a; x) acts 
as a deforming function. 
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3 General procedure 



The first step in our procedure consists in taking for V c g any known SI potential under 
parameter translation. This means that the initial potential V will be determined by 
inverting (3) as V = V cS — V, so that its parameters a will be a = (6, a, £). 

We then consider H as the first member H = H of a hierarchy of Hamiltonians 

i 

Hi = A + (a,\)A-(a,Xi) + Y, e i> « = 0,1,2,..., (5) 

3=0 

where the first-order operators 

A ± (cx, \) = TVf(^,x)^y/f{a;x) + W(X t ; x) (6) 

satisfy a deformed SI condition 

A-(a,\i)A + (a, \) = A + (a, X i+1 )A-(a, X i+1 ) + e i+1 (7) 

for i = 0, 1,2, Here q and X i: i — 0, 1, 2, . . . , are some constants. Solving Eq. (7) 

means that it is possible to find a superpotential W(X;x), a deforming function f(a;x) 
and some constants Àj, e^, i — 0, 1, 2, . . . , with A = A, such that 

Ks(b; x) = W 2 (X; x) - /(a; x) W(A; x) + e (8) 

and 

^(A^ + ZK^W^A^) 

= ^ 2 (A m ;x)-/(a; a ;)Vy'(A m ; a ;) + e m , 2 = 0,1,2,.... (9) 

As a consequence, the (deformed) SUSY partner H\ of H will be characterized by a potential 
V c ff,i(6, a, A; x) = V c g(b; x)+2f(a; x)W'(X; x). Note that in Eqs. (8) and (9), the additional 
terms with respect to the undeformed case are proportional to gW . 

To find a solution for W, f, e« and Àj, our strategy consists in (i) assuming that the 
deformation does not affect the form of W but only brings about a change in its parameters 
A (which will now also depend on a), and (ii) choosing g(a;x) in such a way that the 
function g(oc;x)W'(X;x) contains the same kind of terms as those already present in the 
undeformed case, i.e., W 2 (X;x) and W'(X;x). 

Once we have found a solution to Eqs. (8) and (9), the bound-state energy spectrum 
and corresponding wavefunctions can be found as in conventional SUSYQM. They may be 
written as 

n 

E n (a,X) = J2^ (10) 

i=0 

and 

iíï,^ 1 , x x ( f x W(X n ,x) ^\ 

ïp n (a, X;x) oc = (p n (a, A; x) exp - / — —dx), (11) 

y/f(at;x) \ J J{oí;x) j 
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where <^ n (o;, A; x) fulfils the equation 

<p„+i(a, A; x) = -/(a; :c)<pí,(a, À x ; x) 

+ [W(X n+1 ; x) + W(A; x)]<p n (a, A i; x) (12) 

with <£o(cí, A; x) = 1. They remain however formal solutions till one has checked that they 
satisfy appropriate physical conditions. In the PDM or deformed case, such conditions are 
twofold: (i) square integrability on the interval of definition of V c g, as in the conventional 
case, and (ii) Hermiticity of ir or, equivalently, of H in the corresponding Hilbert space. 
The latter imposes that \ip\ 2 f = ip 2 /VM vanishes at the end points of the interval. This 
extra condition may have some relevant effects whenever the PDM vanishes there. 



4 Classes of superpotentials and corresponding 
g{(X',x) 

One can show that all the known potentials which are SI under parameter translation fall 
into three classes. For all of them, we have determined the general form of the deforming 
function allowing Eqs. (8) and (9) to remain solvable in accordance with the general strategy 
reviewed in Sec. 3. 

In terms of some parameter-independent function <f>(x) and of two parameters À, /i 
making up the set A, the results are given by 

Class 1: W(X; x) = À0(x) + /i, 

(j)'(x) = A(f) 2 (x) + B(f)(x) + C, 

_ A'(cx)cj ) 2 (x) + B'(ex)<P(x) + C'(ex) 
9{ ' ' A(f 2 {x) + B<f>{x) + C ' 



Class 2: W(X; x) = \</>(x) 



JL 



0(X)' 

<t>'{x) = A(f) 2 (x) + B, 
9 (c*;x) = A '^ 2 ^ + B '^ 



(13) 



A(f) 2 (x) + B 
Class 3: W(X;x)= MM±i^ 



^A(t) 2 {x) TW_ 
4>'{x) = [C(j)(x) + D} v /A<j) 2 (x) + B, 

a(a . r) - C'(<*)<Kx) + D'( a ) 
g{OL,x } - C( f)(x) + D 

Here A, B, C, D and A'(a), B'(ct), C'(oc), D' (ex.) are some numerical and a-dependent 
constants, respectively. 

In all three cases, the integral on the right-hand side of Eq. (11) can be explicitly carried 
out by simple integration techniques. Furthermore, by changes of variable and of func- 
tion, ip n (cx, X;x) in the same equation can be transformed into an nth-degree polynomial 
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P n (a, A; y) as follows: 

Class 1: <p n (a, A; x) = P n (a, X;y), y = <j)(x), 

Class 2: ip n (a, X, x) = y~ n/2 P n (a, X, y), y = ( f ) ~ 2 (x), 

Class 3: v?„(a, A; x) = (Ay 2 + B)~ n/2 P n (a, X, y), y = </>(x). 



(14) 



Such polynomials are related to deformed classical orthogonal polynomials and satisfy the 
equations: 

Class 1: P n+1 (a,X;y) = —{[A + A'(a)}y 2 

+ [£? + B\a)]y + C + C'(a)}P n (a, A i; y) 

+ [(A„+i + \)y + Hn+\ + lA p n{ot, Xi,y), 
Class 2: P n+1 (a, X; y) = 2y{A + A'{a) + [B + B'(a)]y}P n {a, A l5 y) 

+{A„+i + A - n[A + A'(a)] + [/x n+ i + // - n(B + B'(a))]y} (15) 

xP n (oc,Xi,y), 
Class 3: P n+ i(a, A; y) = {[C + C"(a)]y + D + £>'(<*)} 

x -(Ay 2 + B)P n (a, Xi,y)+ nAyP n (a, X t ; y) 
+ [(A„+i + \)y + Hn+\ + lA p n{oí, Xi,y). 

5 Some simple examples 

5.1 Particle in a box and trigonometric Pòschl-Teller potential 

Let us consider the superpotential 

W(\;x) = Atanx, ~ < x < (16) 



In the undeformed case, for X = A > 1, it gives rise to the trigonometric Pòschl-Teller 
potential [9] 

14ff(A; x) = A(A - 1) sec 2 x, (17) 

whose bound-state energies and wavefunctions [10] are given by E n = (A + n) 2 and ip n (x) oc 

(cosx) A CÍ" 4 ' ) (sinx), n — 0, 1, 2, The particle-in-a-box problem being the limiting case 

of Pòschl-Teller for A — > 1 corresponds to 



K ff (x) = 



if -f < x < f 
oo if X = ±| 



(18) 



The superpotential (16) belongs to class 1 with (f>(x) = tanx and fi = 0. On choosing 
A = C = 1, B = 0, A'(a) = a, and B'(a) = C'(a) = 0, we get 



g(a; x) — a sin 2 x, 



(19) 
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which for — 1 < a ^ leads to a positive-definite deforming function f(a;x). 

In the particle-in-a-box problem, one easily finds that Eqs. (8) and (9) are fulfilled 

provided Aj = (i + 1)(1 + a) and e; = (2i + 1)(1 + a), i — 0, 1, 2, The corresponding 

energy eigenvalues and wavefunctions are given by 

E n (a, A) = (l + «)(n+l) 2 , (20) 
ícos x) n ~'~^ 

?/>„(«, A;x) oc — r^— — — P n (a, A; tanx), (21) 

(1 + a sm xp n+z » z 

respectively. Here P n (a,X;y) satisfies the equation 

P n+1 (a,\;y) = -[1 + (1 + a)y 2 ]P n (a, A x ; y) 

+ (n + 3)(l + a)yP n (a,X 1 ;y), (22) 

where a dot denotes derivative with respect to y. For any n — 0, 1, 2, . . . , the wavefunc- 
tions (21) are square integrable and ensure the Hermiticity of ir. Hence in the presence 
of deformation (19), the particle-in-a-box problem still has an infinite number of bound 
states making up a quadratic spectrum. As can be checked, Eqs. (20) and (21) go over to 
the undeformed energies and wavefunctions since P n (a, A; tanx) becomes proportional to 
sec™ x Cn^ (sin x) for a — > 0. 

For the trigo nometric Pòschl-Teller potential (17), the results are similar although more 
complicated. The energies and associated wavefunctions are then given by 

E n {a,X) = (A + n) 2 - a{X - n 2 ) 

= [|(A + l)+n] 2 + an(n+l)-ia 2 , (23) 
ÍJ n (a,X;x) oc (cosx)ï^ +n (l + a sin 2 a;)-Kï1^ +n+1 ) 

x P n (a, A; tanx), (24) 

where A = |(1 + a + A), A = a/(1 + a) 2 + 4A(A — 1) and P n (a, A; y) satisfies the equation 

P n+1 (a,X;y) = -[1 + (1 + a)y 2 ]P n (a, A x ; y) 

+ [2A + (n + 1)(1 + a)]yP n (a, A i; y) (25) 

with Ai = A + 1 + a. All functions ip n (a, X; x), n = 0, 1, 2, . . . , are physically acceptable 
as bound-state wavefunctions. 

The starting potential in the PDM Schròdinger equation (1) can be written as V = 
V c g — V, where V e g is given by (17) or (18), while 

V(a, p, a;x) = — (p + a)a 2 cos 2 Ix + pa(2 + ct) cos 2x + aa 2 . (26) 

5.2 Free particle and hyperbolic Pòschl-Teller potential 

In the undeformed case, the hyperbolic counterpart of the superpotential (16), namely 

W(X; x) = X tanhx, — oo < x < oo, (27) 
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corresponds for À = A > to the hyperbolic Pòschl-Teller potential [11] 

V cS (A; x) = —A(A + 1) sech 2 x, (28) 

whose n max + 1 bound-state energies and wavefunctions are given by E n = — (A — n) 2 and 

ip n (x) oc (sech x) "Cn 2 (tanhx), n — , 1, . . . , n max , with A — 1 < n max < A For 
A — > 0, we get the free-particle problem as a limiting case. 

The superpotential (27) also belongs to class 1 with (f>(x) = tanhx and fi = 0. On 
choosing this time A = — C = —1, 5 = 0, A'(ct) = a, and B'(a) = C'(ct) = 0, we get 

g(a; x) = a sinh 2 x, (29) 

which for < a < 1 leads to a positive-definite deforming function f(a; x). 

In the case of the hyperbolic Pòschl-Teller potential and g(a;x) given in (29), the 
energies and associated wavefunctions can be written as 

E n (a, X) = -(A - n) 2 - a(X + n 2 ) 

= -[i(A-l)-n] 2 + an(n + l) + \a 2 , (30) 

ip n (a,X;x) oc (sech x) ï^~ n (l + a sinh 2 x)^^' n ^ 

x P n (a, X; tanhx) (31) 

where X — \(a— 1 + A),A = y/JT— a) 2 + 4A(A + 1) and P n (a, À; y) satisfies an equation 
similar to (25). However, it turns out that although for any n = 0, 1, 2, . . . , ip n (a, À; x) is 
square integrable on the real line, it does not satisfy the condition \ijj\ 2 f — ■> at the bound- 
aries x — > ±oo. Hence, with a deformed function corresponding to (29), the hyperbolic 
Pòschl-Teller potential has no bound state. 

This result can be extended to the free-particle problem. It contrasts with what was 
obtained in [5] in another context and illustrates the strong dependence of the bound-state 
spectrum on the mass environment. 

6 Results and comments 

We have used the procedure illustrated in Sect. 5 to obtain a deforming function satisfying 
Eqs. (8) and (9), as well as the resulting bound-state energies and wavefunctions, for all 
the SI potentials contained in Table 4.1 of [12]. Below we list the results obtained for g, 
E n , and V. 

Shifted oscillator: 

V cS = r [x--j , 

g = ax 2 + 2(3x, a > (3 2 > 0, 
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[{2n + 1)A + (2n 2 + 2n + l)a](5 - bu x 2 
_____ 

A = Veo 2 + a 2 , n = 0,1,2,..., 
V = 2(p + 2a)ax(ax + 2(3) + 2pa + Aa(3 2 . 
Three-dimensional oscillator: 



, , , 1(1 ■ L) 
V cS = -cu x + — — ^— L , < x < oo, 

x 2 

„2 



í X' 

g = ax 2 , a > 0, 

2(n + / + l)(2n+l) + ^ 



a 



Coulomb: 



E„ — 



E n = A (2n + l + + 

A = Veo 2 + a 2 , n = 0,l,2,.. 
V = 2(p + 2a)_V + 2p_. 



^ = __ + _+_!, 0<x<oo, 
re ar 

g = «ï, et > 0, 
e 2 -a[n 2 + (/+l)(2n + l)] x 2 



2(n + / + 1) 

n = 0, 1, . . . , n max , where n max = largest integer such that 
n 2 + (/ + l)(2n+l)<Çifa<^, 



Morse: 



V = aa 2 . 



V eS = B 2 e~ 2x - B(2A+l)e~ x , A,B>0, 
g = ae~ x , a > 0, 

1 / 25(2,4 + 1) - [(2n + 1)A + (2n 2 + 2n + l)a] x 2 
n " ~4 V A + (2n + l)a 



A = V45 2 + a 2 , n = 0, 1, ... , n max , where n max = 
largest integer smaller than A and such that a < a max (n max ) with 
AA(A + l)B 



ttmax(O) 



2A + 1 ' 
B(2A+ l)(2n 2 + 2n+ 1) 



2n 2 (n + l) 2 
B(2n + 1)[(2A + l) 2 + An 2 (n + l) 2 ] 1 / 2 



2n 2 (n+l) 2 ' U 1 ' 2 '···' 



V = (p + a)a 2 e 2x + pae~' J 



Eckart: 



V eS = A(A - 1) csch 2 x -2B coth x, A > -, B > A 2 , 
< x < oo, 

g = ae~ x sinhx, -2 < a ^ 0, 



Er, 



-{A + nf- 



B - \a[{2n + l)A + n 2 } 



A + n 
-a[(2n+ l)A + n 2 ], 
n = 0,1,2,... if a = -2, 

n = 0, 1, . . . , n max if a > — 2, where n max = largest integer 
such that (A + n) 2 < 2B+ °ff~ 1} , 



K = (p + o-)a 2 e~ 4:r - pa(2 + a)e 



-2x 



Scarfí: 



V eff = [B + A — A) sec x - B{2A - 1) tanx sec x, 

< B < A-l, -- < x < -, 

2 ~ - 2 

g = ct sinx, < \a\ < 1, 



-En 



-I(2n+l + A + + A_) 2 + a(n + i ) (A+ - A_) 
- a 2 (n 2 + n + ± ) 



A ± = ^(l T a) 2 + (A±5)(i±5-l), n = 0,l,2,. 



V = — (p + er)» sin x — pa sinx + aa . 



Rosen-Morse I: 



V cS = A(A - 1) esc 2 x + 2B cot x, A > -, < x < vr, 



g = sin x(ct cos x + /3 sinx), ^ < a/1 + /3 > — 1, 
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E n = (A + n) 2 -í 2 A + J -) +(3[(2n+l)A + n 2 }, 

n = 0,l,2,..., 

F = (p + a) 

+ p(2 + /3)(-asin2a; + /3cos2a;) + (-p + a)Ua 2 + (3 2 ). 

Three potentials considered in [12] are missing from the list: Scarf II because no non- 
trivial vàlues of the parameters may ensure positive definiteness of f(ac;x), Rosen-Morse 
II and generalized Pòschl-Teller because they do not have any bound state in the deformed 
case. 

For the remaining potentials, strikingly distinct influences of deformation or mass pa- 
rameters on bound-state energy spectra are observed. In some cases (shifted oscillator, 
three-dimensional oscillator, Scarf I and Rosen-Morse I), the infinite number of bound 
states of conventional quantum mechanics remains infinite after the onset of deformation. 
Similarly, for Morse potential and for Eckart potential with a / -2, one keeps a finite 
number of bound states. For the Coulomb potential, however, the infinite number of bound 
states is converted into a finite one, while for Eckart potential with a = —2, the finite 
number of bound states becomes infinite. It is also remarkable that, whenever finite, the 
bound-state number becomes dependent on the deforming parameter. 

For the potentials V to be used in the PDM Schròdinger equation (1), we get either the 
same shape as V c ff (shifted oscillator, three-dimensional oscillator, Coulomb and Morse) or a 
different shape (remaining potentials). In the first case, the mass and ambiguity parameters 
only lead to a renormalization of the potential parameters and/or an energy shift. 

7 Conclusion 

In this communication, we have shown how to generate new exactly solvable PDM (resp. de- 
formed) Schròdinger equations with a bound-state spectrum different from that of the corre- 
sponding constant-mass (resp. undeformed) Schròdinger equations and we have illustrated 
our method by several examples. In addition, we have demonstrated the importance of 
the Hermiticity condition on the deformed momentum operator for the existence of bound 
states. 
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- (a 2 — j3 2 ) cos 4x + af3 sin 4x 
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